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Abstract 

Here, in order to study stochastic resonance (SR) in a double- well potential when the noise source 
has a spectral density of the form f~ K with varying k, we have extended a procedure, introduced 
by Kaulakys et al (Phys. Rev. E 70, 020101 (2004)). In order to have an analytical understanding 
of the results, we have obtained an effective Markovian approximation, that allows us to make 
a systematic study of the effect of such kind of noises on the SR phenomenon. The comparison 
of numerical and analytical results shows an excellent qualitative agreement indicating that the 
effective Markovian approximation is able to correctly describe the general trends. 
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I. INTRODUCTION 



Stochastic resonance (SR) is one of the most interesting noise-induced phenomena, that 
arises from the interplay between deterministic and random dynamics in a nonlinear system 
Q]. This phenomenon has been largely studied during more than two decades due to its 
great interest not only from a basic point of view but also for its technological interest and 
biological implications 

HQ. 

Most of those studies have used white or colored noises, with a few exceptions where 
more wide classes of noises were considered. For instance, in |yj SR in systems subject 
to a colored and non Gaussian noise was studied. However, there are other works more 
tightly related with the present one, as indicated by the following examples. In jj], the 
authors studied, through numerical simulations, the behavior of the signal-to-noise ratio 
(SNR) gain in a level crossing detector and a Schmidt trigger, when subject to a colored 
noise composed of a periodic train pulse plus a Gaussian f~ K noise with variable k. Their 
results indicate that the maximum of the SNR is larger for white noise, and moves towards 
ia rg e n 01s e Cities for increasing , In fl expe— evidence was fonnd that noise can 
enhance the homeostatic function in the human blood pressure regulatory system. Related 
with the last work, other experimental evidence was found in |6] that an externally applied 
f^ 1 noise, added to the usual white noise, contributes to sensitizing the baroflex function 
in the human brain. In |7[, and in a model of traffic junction of a main and a side road, 
it was found that the effect of a Gaussian f~ K noise with n > shows an overall traffic 
efficiency enhancement. An enhancement of the SR phenomenon in a FitzHugh-Nagumo 



a 



model submitted to a colored noise with f~ K for > k > 2 was found in iSf. In I2f it was 
experimentally demonstrated that an SR-like effect can be obtained in rat sensory neurones 
with white, f~ l and f~ 2 noises, and that, under some particular conditions, f -1 noise can 
be better than white noise to enhance neuron's response. Related to it, in [l0| it was shown 
that it is possible to enhance the SR effect in a FitzHugh-Nagumo model submitted to a 
colored noise with f~ K , and that the optimal noise variance of SR could be minimized with 
K w 1. 

Motivated by the work of Kaulakys and collaborators , who have introduced a method 
to generate f~ l noises over a wide range of frequencies (see also 12,[l3|), here we discuss 
how to extend such a procedure for positive and negative values of the stochastic variable. 
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We also exploit this procedure to analyze the effect of a noise spectrum of the form f~ K 
with varying k, on the SR phenomenon in a simple double-well potential. In the following 
Section we present the model system to be studied and the procedure to generate the f~ K 
noise. Afterwards we discuss an effective Markovian approximation, and exploit it to study 
SR. Finally we discuss the results and draw some general conclusions. 



II. THE SYSTEM 



A. Stochastic differential equations 

The starting point of our analysis is the following system of stochastic differential equa- 
tions 

x = f{x)+g{x)y{t) (1) 
u(y) D , . . . . . . 

y = — + -v(y)£(t), (2) 

r r 

where x is the coordinate of a particle diffusing in a double well potential U (x) = 
~ f X f(()d( = t — T' su bject to a noise y(t). The second equation corresponds to the 
Langevin equation driving the noise y(t), inspired in the work of Kaulakis et al. [111 ]. In this 
last equation we consider a new potential V(y) = — f y u(s)ds, and a (white) noise that 
enters in a multiplicative form with a function v(y). The last function will be characterized 
by an exponent \i. 

We consider the following form for the function u(y) 

u(y) = ay 3 - (3y 5 + s(y)y 4 (3) 

where s(y) indicates the sign of y (i.e, — 1 if y < and +1 if y > 0). For the function v(y) 
we adopt 

v{y) = |y|" + c, (4) 

where both, the exponent n and the constant c are positive (> 0). 

The above indicated forms ch ange the symmetry of the potential V(y) and, in addition, 



when compared with the work in 



increases the range of the noise variable from [0, +oo) 



to (— oo, +oo), as is shown in Fig. 1. The parameter c allows for the random variable y to 



adopt negative values when c > \y\ 



A' 




FIG. 1: Symmetric potential V{y) as derived from Eq. (J3J). We have used the following values: 
a = 5 x 1CT 4 and (3 = \. 



B. Characteristic of the noise variable y 

The most relevant aspect of the process y is its power spectral density (PSD) with a 1/ f 
frequency behavior. Kaulakis et al have shown that when c = and fi = 5/2, the noise 
y exhibits a 1/f functionality in a wide range of frequencies. For c > 0, but small, this 
property is still valid, as we show in figures 121 and El 

When the exponent [i changes, the PSD behaves as 1/ f k , with k < 1 for fi < 5/2. We will 
use this property to evaluate the mean-first-passage-time (MFPT), and the signal-to- noise 
ratio (SNR). In particular we have used fi = 3/2, yielding k ~ 3/4. 



III. EFFECTIVE MARKOVIAN THEORY 



In order to be able to obtain some analytical results, we resort to an effective Markovian 
app roximation, analogous to the so called unified colored noise approximation (UCNA) 
Il5| |. Exploiting this approach we are able to find an effective Markovian Fokker-Planck 
equation (FPE) for the probability density P(x,t). The procedure is the following (the 
prime will indicate derivation respect to the variable x). 
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FIG. 2: A typical realization for Eq. © using /i=5/2 and c = 1 x 10 -4 . The PSD of this realization 
shows a 1// behavior, see Fig. 03 

A. Adiabatic procedure 



Deriving Eq. respect to the time we have 

x = f'(x)x + g'(x)xy + g(x)y. 

Now, assuming an adiabatic behavior, we eliminate x, and using Eq. (0) we obtain 

~x - f(x) 



~ f'{x)x + g'(x)x 



+ g{x) 



u(Z(x)) | Dv(Z(x))gt) 



T 



T 



(5) 



(6) 



where we have defined Z[x) = Z (x) + Zi(x), whith Z (x) = — ^fej and Z x {x) = Now, 
we use the following approximation 



u(Z(x)) Pd u(Z (x)) + u'(Z (x)) Z ± (x) 



(7) 



and similarly for v 



v(Z)^v(Z )+v'(Z )Z 1 (x) 



Adopting now g(x) = 1, that implies Zq = —f(x), we have 

o = f(x)x + ^Mlx + D^lew + o(±m- 

T T 



(8) 



(9) 
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FIG. 3: PSD for the variable y, as indicated in Eq. (J2J. We used the same values of parameters as 
in Fig. |2J The white line corresponds to a linear fitting, resulting in a slope k = —1.004 ± 0.005. 

With the above indicated results, the effective equation for the process x adopts the following 
form 



x 



u(Z {x)) + Dv(Z {x))£(t) 



A 1 (x) + B 1 (x)C(t), 



(10) 



rf'(x) + u'(Z (x)) 

and, due to the polynomial character of the function h(y), we can write the following limit 
for A(x) and B(x) when r — > 

Zq(x) 



Ai(x) 
B 1 {x) 



A{x) 
Dv{Z {x)) 



B(x). 



(11) 
(12) 



u'(Z (x)) 

Finally, using the above indicated approximations, the stochastic differential equation for 
the process x reads 



x = A(x) + B(x)£(t). 



(13) 
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B. Fokker-Planck equation 

The FPE associated with the Langevin equation, Eq. JT3J), is (using the Ito prescription 

^P(x,t) = -JL[A(x)P{x,t)) + l^- 2 [B\x)P(x,t)]. (14) 

As is well known, the stationary distribution of this FPE is given by 

N 

P{x) = ^yexp {-*(*)}, ( 15 ) 
where N is the normalization factor, and 

*{x) = 2 f ^±d(. (16) 

C. Mean-First-Passage-Time and SR 

The indicated FPE and its associated stationary distribution allow us to obtain the mean- 
first-passage-time (MFPT) through a Kramers-like approximation. Using known expressions 



we obtain for the MFPT 1 1( 



\ f x ° dy f y dz^(z) 

™ = 2 J. m L-ser- (17) 



vp(i) = exp I 2 f I CIS) 



where 



In order to study SR, as usual we introduce an external signal in the form of a term 
rocking the double well potential: U(x) = U (x) + S(t), with S(t) = S sxn(ujt) (in what 
follows we adopt u> = 1.33 x 10 -5 ). Exploiting the so called "two-state approximation" jlj|, 
we define the SNR as the ratio of the strength of the output signal and the broadband noise 
output evaluated at the signal frequency u, obtaining 



utpi 

a 



SNR <x{if} s . , (19) 
where the derivative of the T in the above expression, as indicated, is evaluated at S = 0. 
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IV. RESULTS AND CONCLUSION 



We have done extensive numerical simulations of the full set of Eqs. (jl!2|) in order to 
obtain the SNR. The results are shown in Fig. 4. Also, in Fig. 5 we show the SNR computed 
using the effective Markovian theory, obtained through Eqs. (fT7|) and (|19J) . In order to be 
able to compare the results we have normalized the curves. Also, in order to have a well 
defined variance of the noise process, in all the simulation we have obtained numerically the 
variance of y(t), as described by the Eq. (J2J). 
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FIG. 4: SNR obtained when simulating the full set of Eq. (11121) . Here a corresponds to the noise 
intensity defined through the distribution width, as indicated in the text. Squares and circles 
corresponds for fi = 3/2 and fi = 5/2 respectively. The lines are for guiding the eye only. 

From the comparison of both figures it is apparent that the results obtained using the 
effective Markovian theory are in very good (qualitative) agreement with those from simu- 
lations. This is in accord with previous results obtained for different systems jsl ^| . We 
can conclude that such kind of UCNAdike approximation offers an adequate framework to 
obtain effective Markovian approximations in a very wider class of systems than the one to 
which was originally applied 

The above indicated results are in complete agreement with those of |4j. That is: the 
maximum of the SNR is larger for white noise, and it moves towards large noise intensities 
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FIG. 5: SNR obtained using Eq. (|19[). as derived from the two state theory. Continue and dashed 
line correspond to (i = 3/2 and fi = 5/2 respectively. 

for increasing k. In order to gain some physical insight about this behavior it is worth to 
remark that the function defined by Eq. ()18|) is directly related to an effective potential 
within the approximation we used 



The behavior of such a potential reveals what are the consequences of changing the exponent 
/i: when /i = 5/2 (i.e. the PSD is 1//), the effective potential shows a well defined well; but 
as /i decreases the well is less defined as shown in Fig. 6. Such a behavior of the effective 
potential explains why the SNR increases when /i decreases. The general theory shows that 
the SNR increase is proportional to the Kramers rate rjc, that is given by 



where AV e // is the high of the barrier in the effective potential separating the attractors. 
Hence, the reduction of the SNR with increasing k (or decreasing fi) could be directly related 
with the marked reduction of the barrier separating the attractors in the effective potential 
picture. This also explains the reason for the shift of the SNR maximum towards larger 
values of the noise intensity. 




(20) 




(21) 
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FIG. 6: Effective potential, Eq l|18[). showing the different behavior of the wells for different values 
of the exponent [i: continuous line \x = 5/2, dotted line fi = 3/2. Due to the symmetry, we only 
show positive values of x. 

In conclusion, we want to remark that this is a first step towards an analytical under- 
standing of two very important, connected, and ubiquitous, aspects in natural processes. 
Those are: the l/f k behavior of noises' PSD, and its role in signal detection via the SR 
mechanism. The physical picture provided by the indicated effective Markovian approx- 
imation, offers an adequate framework to analyze and understand the main qualitative 
trends of such phenomenon. Even more, we expect to apply the same scheme to other noise 
induced phenomena when subject to l/f k noises. This will be the subject of further work. 
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